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Abstract 

We explore the viability of a boson dark matter candidate with an asymmetry between the 
number densities of particles and antiparticles. A simple thermal field theory analysis confirms 
that, under certain general conditions, this component would develop a Bose-Einstein condensate 
in the early universe that, for appropriate model parameters, could survive the ensuing cosmological 
evolution until now. The condensation of a dark matter component in equilibrium with the thermal 
plasma is a relativistic process, hence the amount of matter dictated by the charge asymmetry is 
complemented by a hot relic density frozen out at the time of decoupling. Contrary to the case 
of ordinary WIMPs, dark matter particles in a condensate must be lighter than a few tens of eV 
so that the density from thermal relics is not too large. Big-Bang nucleosynthesis constrains the 
temperature of decoupling to the scale of the QCD phase transition or above. This requires large 
dark matter-to-photon ratios and very weak interactions with standard model particles. 

PACS numbers: 98.80.-k, 98.80.Jk, 95.35.+d, 95.30.Sf, 03.50.-z 
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I. INTRODUCTION 


The idea of a weakly interacting massive particle (WIMP) has guided the research in 
dark matter (DM) over the last three decades [TH3]. In the standard cosmological scenario 
WIMPs are produced thermally in the early universe, with an abundance of particles today 
fixed at the time of thermal (chemical) decoupling. The energy density in WIMPs is mainly 
determined by their self-annihilation rates into the standard model sector. For particles with 
masses and couplings at the electroweak scale the energy density in thermal relics roughly 
coincides with the known abundance of DM. Particles with these properties appear in e.g. 
supersymmetric extensions to the standard model, and this coincidence is usually known as 
the “WIMP miracle.” 

However, there are alternatives that deserve to be considered seriously. We know that 
the amount of baryons in the universe is not determined at thermal decoupling, but by an 
asymmetry between the number densities of particles and antiparticles. Something simi¬ 
lar could have also happened to DM [3], If in addition to the existence of a nonvanishing 
conserved charge, DM is described in terms of a boson particle that was in thermal equi¬ 
librium with the constituents of the standard model in the early universe, it is natural to 
think that the zero-mode could have developed a nonvanishing expectation value, leading 
to the appearance of a cosmological Bose-Einstein condensate (BEC) [•") 7|. We call this 
scenario asymmetric condensed dark matter. Other proposals with a similar spirit have been 
previously considered under many different names [SfKS]. See also Refs. |2ZH3U for similar 
ideas in the context of the QCD axion, where unlike the case considered here DM never 
thermalizes with the standard model particles in the early universe. 

In this paper we identify four conditions that, together, guarantee the appearance of a 
cosmological thermal condensate of particles that could match the totality of the cold dark 
matter (CDM) today: 

1. DM is described in terms of a boson field with a conserved current, 

2. There is a nonvanishing total boson charge with the right amount in the universe, 

3. The boson was in thermal (kinetic and chemical) equilibrium with the particles in the 
standard model, 


2 


4. Thermal decoupling took place at an energy scale below the critical temperature of 
condensation. 

Under these assumptions we conclude that a thermal BEC can only emerge from a rela¬ 
tivistic process, and that the particles composing the condensate should be much lighter 
than standard WIMPs. This requires DM-to-photon ratios that are at least eight orders of 
magnitude larger than in the case of baryons, a point that deserves further attention. It is 
this large asymmetry what makes it possible that light particles are nonrelativistic in spite 
of being in thermal equilibrium in the early universe. The interaction of the DM bosons 
with the standard model particles must be very weak, even more than the weak nuclear 
force, although strong enough to reach early thermal equilibrium. We speculate about the 
possibility that the gravitational interaction alone could bring the DM bosons into thermal 
equilibrium with standard model particles, but the naive analysis we present in this paper 
deserves a further study in terms of the Boltzmann equation. Figure [l] in Section III 


sum¬ 


marizes the main features of this proposal, that remain valid if the condensate contributes 
only to a partial fraction of the total CDM in the universe. 

Note that the four conditions above lead to a macroscopic excitation of the zero-mode 
that in the present universe is not precisely a thermal BEC by the standard definition, i.e. 
it has not been in thermal equilibrium since DM decoupling. We will nonetheless refer to 
this state as the condensate in this paper. We should emphasize that the description of a 
thermal condensate we present here is only valid as long as the expansion rate of the universe 
is lower than the characteristic frequency of oscillation of the zero-mode. Furthermore, in 
order to not affect the abundance of light elements predicted in the standard cosmological 
scenario the DM particles should be already decoupled from the thermal plasma before the 
outset of Big-Bang nucleosynthesis (BBN), and then we cannot extrapolate the results in 
this work to candidates lighter than around m ~ 1CD 14 eV. Further constraints on the tem¬ 
perature of decoupling coming from the effective number of extra neutrino species increase 
this value to 10~ 12 eV. Particles with a mass below this number require a more detailed 
analysis that we leave for another paper. We also stress that our motivation here is purely 
phenomenological, and we do not develop any specific boson candidate (a possible particle 
realization is discussed in e.g. Ref. [32], see also m- Rather, we assume an asymmetry 
and an early stage of thermal equilibrium, and explore the consequences of demanding a 
late-time coherent behavior on the model. 
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The paper is organized as follows. In Section [TT] we review some basics on quantum fields 
in curved spaces, and frame the idea of a cosmological thermal condensate of scalar particles 


within our current understanding of physical theory. Then in Section III we use some basic 
cosmological observations to constraint the parameters of the boson as a model of DM. 


Finally, we conclude in Section IV with a brief discussion of the main results of this paper. 


II. COSMOLOGICAL GRAVITATING BOSONS IN THERMAL EQUILIBRIUM 


Let us first consider the general problem of a cosmological gravitating condensate. In this 
paper we assume that the spacetime metric is a classical object, but matter is described in 
the language of quantum field theory. This picture is expected to be appropriate as long as 
we are only interested in describing physical processes taking place well outside the Planck 
regime, with the classical spacetime and the quantum matter fields connected through the 
semiclassical Einstein field equations, G^ u = 8ttG(T^ u ). Here is the Einstein tensor, G 
denotes Newton’s gravitational constant, and (T/j,) is the expectation value of the energy- 
momentum tensor operator, discussed further below. 

In this section we consider the case of a cosmological gas of bosons in thermal equilibrium. 
At this point we do not need to specify how this thermal regime could have been achieved 
(and even supported) in practice. The process of DM decoupling, essential for any sensible 


condensate model, is discussed next in Section III 


A. Quantum field theory in curved spacetimes: a minimal review 

We restrict our attention to the case of a massive, complex scalar field with an internal 
1/(1) global symmetry and no self-interactions. Associated with the global symmetry there 
is a conserved current, and then (crucially) the possibility to have a non vanishing conserved 
charge. The generalization to more involved fields and/or symmetries is possible, but we do 
not consider it here. Classically, the scalar field ip(x) satisfies the Klein-Gordon equation, 

(□ — m 2 )(p(x) = 0 , (1) 

with □ = the d’Alembert operator in four dimensions and x a generic point in 

the spacetime manifold. A minimal coupling between the scalar field and the gravitational 
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interaction has been considered here. We are using units with c — h— 1, where the scalar 
held <p{x) and the mass parameter m have dimensions of energy. 

Given two solutions p\(x) and <p 2 (x) to the classical equation of motion, Eq. (JTj) , we can 
define the symplectic (or Klein-Gordon) scalar product 

(>i, ¥> 2 ) = ~i J [<P 1 (^ 2 ) ~ (d^i)(p* 2 \ . (2) 

Here is a timelike, future-directed, normalized four-vector orthogonal to the three- 
dimensional Cauchy hypersurface E, and dT, = ^yd 3 x is its volume element. The function 
<p 2 (x) in Eq. (J2]) denotes the complex conjugate of <p 2 (x). By construction the Klein-Gordon 
scalar product does not depend on the Cauchy hypersurface E, and thus is preserved with 
time evolution. 

At the quantum level the state of the system is described in terms of a vector in a 
Hilbert space, with the held, <p(x), and its conjugate momentum, tt(x) = (n fJj d^ l p*(x)), 
promoted to operators acting on the elements of this space [MljdT] . In this paper we follow 
the canonical approach and impose the standard equal-time commutation relations 

[<p(t,x),n(t,y)] = i5(x-y), [(pit, x), ip(t, yj] = [n(t, x), tt(£, y)\ = 0. (3) 

We are working in the Heisenberg picture, where the state vectors remain hxed and only 
the held operators evolve in time. As usual we can decompose <p{x) in terms of the time- 
independent creation and annihilation operators, 

fa) = , ( 4 ) 

i 

with Ui{x) a complete set of normal modes, each solving Eq. ([!]), that are orthonormal with 
respect to the symplectic product: 

(ui, Uj) = S i:j , {u*,u*) = —Sij , (ui, u*) = 0 . (5) 

Here the letters i and j are just labels for the mode-functions, and not spacetime indexes. 
The normal modes Ui(x) in Eq. (Jdj) have dimensions of energy, with the creation, aj, b\, and 
annihilation, di, bi, operators dimensionless. All these conventions guarantee the standard 
commutation relations [a*, a]] = [&?:,&]] = S t j , with all other combinations zero. Note that 
in the continuum limit the sum in Eq. Q should be replaced by an integral, and then the 
Kronecker delta in Eq. ([5]) by a Dirac delta function, but the construction is similar. 
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The vacuum is defined to be the state for which a,|0) = 6*10) = 0 for all i. The Hilbert 
space can then be constructed (a la Fock) by successive applications of creation operators 
on the vacuum state. Incidentally, the elements of this construction, |7V“,..., Nf,.. .), are 
the eigenstates of the particle (iVf = d\di) and antiparticle (iYf = b\bi) number operators. 

If the state of the matter fields is described by a vector |-0) in Hilbert space (i.e. a 
pure state), then (0(x)) = ('i/\0(x)\'i/>), with 0(x) denoting a generic observable in the 
theory. However, we wish to describe boson particles at finite temperature, k B T = 1/(3, 
where k B = 1 is the Boltzmann constant. In this latter case we should replace the standard 
expectation value (O(x)) above by [38, 39] 

Tr [e-^ 6 ~^d( X )\ 


{0(x))e = 


( 6 ) 


Tr[e~P(H-vQ)] 

Here Tr[e _|S ^ _/i ^], the partition function, represents the trace of the density matrix in the 
grand canonical ensemble, with H the Hamiltonian operator for the scalar field and p the 
chemical potential associated to the conserved charge Q = —i J[(fk — jT^fl]d 3 x. This charge 
essentially represents the difference between the number of particles and antiparticles in the 
configuration. 


* * * 

Next we derive a number of cosmologically-relevant expressions for this model. Note that 
many of these equations would also apply if there are no conserved charges in the dark sector, 
H — 0, in which case it is not possible to develop a thermal condensate of DM particles. For 
a matter component described in terms of a real field particles and antiparticles coincide, 
and we cannot define properly a conserved charge operator Q in the system. However, if this 
field were not coupled to the standard model of particle physics (or if the interactions with 
this sector were weak enough), it could be possible to reach a diluted state of nonrelativistic 
kinetic equilibrium where, at the effective level, the number of particles were conserved. 
Note however that this state would not be in chemical equilibrium with the cosmological 
thermal plasma, as we assume in this paper. 

It has been recently proposed that this could happen for the QCD axions [23128], where 
the particles are generated initially in a low momentum state through a misalignment of the 
vacuum in the early universe |4(TI 144] , although there is still some debate in the literature [29- 
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[3lj| . In any case, the temperature in axions would not coincide with the temperature of the 
standard model particles, and then the results in this paper do not necessarily apply there. 


B. Cosmological considerations 


In this section we apply a simple analysis in field theory to describe an asymmetric 
thermal gas of bosons in an expanding universe. We concentrate on the appearance of a 
cosmological condensate in the limiting case of a relativistic regime. As we will find next in 
Section |III| the condensation of a DM component in equilibrium with the thermal plasma 
is necessarily a relativistic process, and then a relativistic analysis is mandatory for the 
purposes of this paper. 

In absence of gravity, the emergence of a thermal condensate has been extensively con¬ 
sidered in the literature pOTIBiZ] . Let us now focus on the case in which the boson particles 
live on a spatially flat, homogeneous and isotropic universe. Under this assumption the 
spacetime metric can be written in the form 

ds 2 = — dt 2 + a 2 (t)(dx 2 + dy 2 + dz 2 ). (7) 


Here t is the cosmological time, a(t ) the scale factor, and ( x , y , z) a set of comoving spatial 
coordinates. Without any loss of generality we can choose a = 1 today, when comoving and 
physical quantities coincide. To proceed, we require a set of normal modes adapted to the 
cosmological background. Because there is no timelike Killing vector field associated to the 
line-element in Eq. (J7]) , there will not exist a set of normal modes in the form of ordinary 
plane-waves, and most of the results in flat spacetime do not formally apply. However, if 
the expansion rate of the universe is not too large, H <C m, we can still write 


M x ) = 


VW> 


a°u k 


exp ( — i / u k dt + ik ■ x 


( 8 ) 


where H = a/a is the Hubble parameter, and co k = a Jm 2 + a~ 2 k 2 the dispersion relation of 
the scalar particle. We are imposing periodic boundary conditions over a box of comoving 
size L, with k = (2ir / L)(n x ,n y ,n z ), V = L 3 , and n x ,n y ,n z = 0,±1,±2,.... Note that the 
wavenumber k labels the different mode-functions in Eq. (|8]), and that the mode-energies 
lO k depend only on k = \k\] this is a consequence of the symmetries associated to the slices 
of constant cosmological time t = const. However, the spacetime background is not static, 
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and then the dispersion relation is time-dependent. Remember that we are interested in 
describing a condensate in an infinite homogeneous and isotropic flat universe, and then we 
can switch to the continuum, = V/(2%r) 3 f d 3 k, if necessary. 

The functions in Eq. (|8| satisfy the Klcin-Gordon equation, Eq. ([!]), and the simplectic 
normalization condition, Eq. ([5]), to the second order in an adiabatic WKB approximation, 
H/oj(k ) <C 1. Then, as long as the the expansion rate of the universe is lower than the 
characteristic frequency of oscillation of the zero-mode, H < m, the expression in Eq. ([8]) 
describes the whole spectrum of mode-functions to the next-to-leading order in the adiabatic 
approximation. In this paper we will be interested only in a description of the condensate 
to the zeroth order in this series expansion, and then we will not keep factors of H/uj{k ) in 
the final expressions. 

We are mainly interested in obtaining cosmological thermal expectation values. First, 
however, we need to clarify some issues related with the vacuum of the theory. Formally, 
the expectation value of a quantity that is quadratic or higher order in field operators di¬ 
verges, even when it is evaluated on the vacuum state. See for instance the Hamiltonian 
as a particular realization of a quadratic operator with ultraviolet divergences. We should 
then use some regularization/renormalization prescription in order to obtain, from these 
quantities, well defined objects when acting on physical states. There are different meth¬ 
ods to renormalize higher order operators, such as proper-time regularization, dimensional 
regularization, zeta-function regularization, point-splitting regularization, or adiabatic reg¬ 
ularization [37]. For the case of a homogeneous and isotropic universe the simplest choice 
is probably adiabatic regularization. This prescription was employed for the first time by 
Parker in [53], and we can use this technique to identify sensible vacuum expectation values 
from ill-defined higher order operators. While important on their own, these vacuum contri¬ 
butions are not really relevant for the problem we want to address in this paper, and we will 
simply omit them in all the expressions below (see for instance Refs. [5T1 [55 J for a renormal¬ 
ized energy-momentum tensor in the context of an expanding homogeneous and isotropic 
universe). In practice, this is essentially equivalent to impose standard normal ordering, e.g. 
: OiO? := a%, : til := b% on higher order operators, and ignore the divergent zero-point 
contributions to the expectation values. With all these assumptions the Hamiltonian and 


( 9 ) 


the scalar charge in Eq. (Jh]) simplify to 

H = '£{N2+Nl) Ut , Q = '£(NI-Nl), 

k k 

where, as mentioned above, we have omitted nondiagonal terms in the expression for the 
Hamiltonian operator that are suppressed by factors of H/oj(k) in the series expansion. For 
completeness, the operator representing the total number of particles (including antiparti¬ 
cles) is by definition given by N = JT(iV“ + Nf). 

Since the Hamiltonian and the scalar charge commute, [H,Q] = 0, we can choose a 
basis of the Hilbert space in which both H and Q are diagonal. To the zeroth order in the 
adiabatic expansion the Hamiltonian and the scalar charge also commute with the number 
operators, [NpH] = [Ni,H] = [N?,Q] = [A r |, Q] = 0, and we can then use e.g. the basis 
obtained using the Fock construction, where all the elements have a def ini te number of 
particles and antiparticles. This choice simplifies the analysis, since only the terms in the 
diagonal of an operator can contribute to the expectation value (0(x))p in this basis. This 
is what happens, for instance, in the case of the energy-momentum tensor, where in practice 
we can just forget the other terms in T )JiV (x) different to (NS + Ni)T}£)(x). Here T$(x) is 
the energy-momentum tensor corresponding to an excitation in the h-mode, 

T$(x) = dpu%d„u\ + dnUjdvUz - g /w (d a utffu% + m 2 u^u^} , (10) 


and we have again neglected those terms suppressed in the adiabatic expansion. 

After some algebra (see for instance Chapter 2 in Ref. [38j for some details in Minkowski 
spacetime), we obtain the thermal expectation values associated to the number density, 
np = (N)p/(a 3 V), charge density, qp = ( Q)p/(a 3 V ), energy density, pp = -(Tff(x))p, and 
pressure, pp = (T™(x))p (no sum in n), of a gas of bosons in a cosmological background, 

1 r°° 

np = « + n h Q ) + J ^ (N% + Nl)k 2 dk , (11a) 

1 C°° 

qp = (nl - nl) + — J ( N a k - N*)k 2 dk , (lib) 

1 r°° 

pp = (no + n b 0 )m + j ^ (N£ + N b )uj k k 2 dk , (11c) 

1 /'°° _ n ~ 2 k 2 

» = 2^l^ + Rb ^ k2 dk- < lld > 

It is easy to see that in this coordinate system the other components of the energy-momentum 
tensor vanish by symmetry considerations. Here N°‘ and N b are the mean occupation number 
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of particles and antiparticles in the different energy levels, respectively, 


N a k = 


K = 


( 12 ) 


exp[/3(u> k -//)] - 1 ’ ' exp \f3{u: k + p)] — 1 ’ 

with n k = N k /(a 3 V) and n b k = N k /(a 3 V). Since particles and antiparticles are in chemical 
equilibrium their corresponding chemical potentials are equal in magnitude but opposite in 
sign. From now on we will omit the subindex /3 and the overbars in order to simplify the 
notation unless this may lead to confusion. 

Some comments are in order here. To proceed, we have taken the L —y oo limit in the size 
of the comoving normalizing box, and approximated sums over energy states by integrals 
over wavenumbers, Y2 k = (E/27T) 3 J The integral terms coincide with the standard 
textbook expressions for the number density, charge density, energy density, and pressure of 
a gas of bosons at finite temperature (see e.g. Section 3.3 in Ref. [TJ for an example in the 
cosmology literature). However, the density of momentum states g(k ) = (V/2n 2 )k 2 vanishes 
at k — 0, and then the zero-mode is not contained in the previous integrals. That is the 
reason for which it has been necessary to write no = ng + Uq, qo = —n b 0 , po = (ng + n b 0 )m. 


and po = 0 (overbars omitted) explicitly in Eqs. (11). As we will show next, these additional 
terms coming from the zero-mode could be relevant (and even dominate) at low temperatures 
if the charge density is different from zero. 

In order to see this note that, since the total charge is conserved in a comoving volume, 


the expression in Eq. (lib) (implicitly) fixes the value of the chemical potential. As usual, 
the condition of a positive number of particles and antiparticles in the different energy levels 
demands /j < m. The mean occupation number is always larger in the zero-mode than 
in any other level, no > n k ^ o, but usually negligible when compared to the total particle 
density, no n. However, if at a temperature different from zero the chemical potential 


approaches the critical value at |/i| —> m, a macroscopic amount of charge qo < \q\ must be 
accommodated in the lowest energy mode, k = 0. A phase transition takes place, signifying 

J. (Note that if we set (5 fy 0, 


the appearance of a BEC 


= m in Eq. (12) the 


expression for no formally diverges, but of course in practice the value of the chemical 
potential is close but not equal to ±m.) 

At this point the field operator (p(x) describing the fundamental degrees of freedom in 
the gas develops a nonvanishing expectation value, fy 0, with the mode of lowest 
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energy described in terms of a coherent pure state in Hilbert space, 


&oW) = e ld \/Vqo\ip) , %^ o |-0) = %|V>) = 0, (13) 

i.e. an eigenstate of the annihilation operator do with eigenvalue e l6 ^/Vqo- Here, as usual, 
we have normalized this state to unity, (fi>\ -0) = 1, the phase 6 is arbitrary, and (without any 
loss of generality) we are assuming a universe with more particles than antiparticles, q > 0, 
i.e. fi G [0, m). We can then neglect the mean occupation number of antiparticles in the 
ground state, n b Q <C ng, and approximate go = '^o — n\ ~ n q . 1 The nonvanishing expectation 
value 

v(x) = {i>\ (k u k( x ) + *>l u *k( x )) W) = \j ( 14 ) 

k 

plays the role of an order parameter for the BEC. Note that this order parameter satisfies 
the Klein-Gordon equation (in the nonrelativistic limit the Gross-Pitaevskii equation for a 
system of non-interacting particles). We can then describe the condensate in terms of a 
classical field theory —classical from the point of view of a held theory, purely quantum 
from a particle interpretation. Introducing the expression for the order parameter into the 
equations for the number density, charge density, energy density, and pressure associated to 
a classical scalar held, we recover the zero-mode contributions n 0 = n g, go — Po — 
and po = 0 previously identified in Eqs. (where we have already neglected the mean 
occupation number of antiparticles in the lowest energy level). 

In practice, the critical temperature of condensation, /3 C , is defined (implicitly) in terms 
of the cosmological charge density, g, through the equation 

g = qth\P = Pel* = m]. (15) 


Here the subindex th makes reference to the fact that only the excited thermal modes, 


k 0, have been included in this expression, qp = g 0 + qth■ The right hand side in Eq. (15) 


1 Note that the results in this paper are not really sensitive to the quantum coherence of the state | ip). 
If the particles were just accommodated in the mode with lowest energy with arbitrary phases (i.e. if 
this mode were described in terms of an eigenstate of the particle number operator with eigenvalue Vqo, 


Ng\i/>) = Vqo\4>)), they would source the same energy momentum-tensor as the state in Eq. (13). However, 


contrary to the case in Eq. (141, the expectation value of the field operator would vanish, (fi(x)) = 0. 
This could have important observational consequences when discussing e.g. the direct detection of this 
matter component, but this is not the subject of this paper. 
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codifies the maximum amount of charge q per unit volume that, at a given temperature 
(3 = {3 C , can be accommodated in the thermal modes once the chemical potential has already 
reached its critical value at n — m. Below this temperature the chemical potential cannot 
grow anymore, and a macroscopic amount of charge qo = Q ~ Qth[/3 > /3 C , P = m] must be 
accommodated in the ground state. 


In general, Eq. (15) should be evaluated numerically. However, this will not be necessary 


for the purposes of this paper. As we will find next in Section HI for the case of a cosmolog¬ 
ical condensate of DM particles thermal equilibrium must be broken during the relativistic, 
high charge density regime, where analytic expressions do exist. 

The relativistic condensation of a gas of bosons in flat spacetime was considered for the 
first time in a seminal paper by Haber and Weldon [50]; see also Refs. [SH [52] for more recent 
works. Now we repeat this analysis for the case of a homogeneous and isotropic universe 


in a regime of slow expansion, H <C m. The functions N£k 2 and N^k 2 in Eqs. (11) are 
peaked around comoving wavenumbers (3a~ 2 k ~ 1, and then only those values of k con¬ 
tribute significantly to the integrals at different times in cosmic evolution. In the relativistic 
regime, (3m -C 1, and to the lowest nonvanishing order in m/(a“ 1 fc), we can approximate 
the dispersion relation in Eq. ([8]) to u>k = a _1 /c. Introducing this expression into Eqs. 0. 
we obtain 


n = («n + n„) + 


1 


7T 2 /3 3 

1 


[Li 3 (e^)+Li 3 (e-^)] , 


3 = K-^o) + ^[Li3(e^)-Li 3 (e-^)] , 

[Li 4 (e^)+Li 4 (e-^)] , 


P — K + n 0 )m + 
1 


7 t 2 (3 a 


V = 


7r 2 /3 4 


[Li 4 (e^)+Li 4 (e-^)] . 


(16a) 

(16b) 

(16c) 

(16d) 


Here Li a (x) denotes the polylogarithm function of order a and argument x. Some limiting 
values that will be of interest soon are given in Appendix [A] From the thermal nonzero- 


modes in Eq. (16b) and the definition of the critical temperature of condensation in Eq. (15), 
we can read 


T r — \ — 


m 


(17a) 


where we have made use of some of the expressions in Appendix [A] in order to arrive to the 
final result. Note that the presence of particles as well as antiparticles is relevant at high 
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temperatures, and then the expression for the critical temperature of condensation should 
be modified with respect to that in standard textbooks on statistical mechanics, see e.g. 
Ref. [05]. The integral terms in Eqs. only counts the excited thermal modes; for the 
configurations below the the critical temperature of condensation a macroscopic amount 
of charge must be accommodated in the lowest energy zero-mode. We should then write 
q = qo + qth[P, P = m], and the charge density in the zero-mode can be evaluated from 


Qo{T) = q 


, T 

1 ~~ T r 


if T <T r . 


(17b) 


As long as g > m 3 we can guarantee a relativistic temperature of condensation, T c m, 


and then the approximations in Eqs. (16) are justified. 

Note that for a gas composed of particles with vanishing mass, the critical temperature of 
condensation grows to infinity. This might seem to suggest that, in the zero-mass limit, any 
charge in the universe should be necessarily accommodated in the ground state. However, 
remember that the presentation in this paper is only valid as long as the expansion rate of 
the universe is not too large when compared to the characteristic time of oscillation of the 
zero-mode [this comes from the choice we made in Eq. Q for the mode-functions that solve 
Eqs. (|TJ) and ([5])], and only if the gas is in thermal equilibrium. In particular, we cannot 
extrapolate the expressions in this section to the massless case. As we will discuss later 


in Section |III[ in order have a successful cosmological picture any sensible DM candidate 
should be already decoupled from the thermal plasma at BBN, and then some caution is 
necessary even for non-zero mass particles; we will give some numbers soon. 

The expressions for the energy density and pressure in the thermal gas are quite illumi¬ 


nating. We can obtain these two quantities expanding Eqs. (16c) and (16d) around (3p • Cl. 
After some algebra we get 


7 r 2 T 4 
~15~ 


p = qo(T)m 


7r 2 T 4 

~ 45 ~ 


P = 


(18) 


Note that there are two different contributions to the expressions in Eq. (18). On the one 
hand that coming from the particles in the condensate, p = qo(T)m, p = 0, that of course 


coincides with the contribution associated to the order parameter in Eq. (14). The transition 


to a condensate is quite sharp: according to the expression in Eq. (17b), if the temperature is 
only one order of magnitude lower than the critical temperature of condensation, then 99% 
of the charge will be already accommodated in the lowest energy level. As a consequence of 


13 













this we can think of qo{T) as a function that steeply transits from zero to the total charge 
density q when the system approaches the critical temperature of condensation from above. 
In addition, and as a consequence of the high temperature in the gas, (3m 1, there is 

relativistic gas of thermal bosons, p = 7r 2 T 4 /15, p = 7r 2 T 4 /45, that appears even in the 
presence of a condensate when the temperature of the universe is lower than the critical 
temperature. This thermal cloud will be essential for the analysis in the next section. 

That the coherent excitation of the zero-mode appears only below a critical temperature 
of condensation might suggest that it is a late-time phenomena, characteristic of an old 
and cold universe. However, the critical temperature of condensation is not constant in 
cosmic history, and changes with the expansion of the universe as T c r^j 1/a 3 / 2 . Then, the 
critical temperature of condensation decreases with the cosmological expansion faster than 
the temperature of a universe in adiabatic expansion, where T ~ 1/a (remember that it 
is the temperature of the CMB photons that determines the temperature of the gas, as 
long as they are in thermal equilibrium due to the self-annihilation of DM particles into the 
standard model sector). Therefore, given an asymmetry between the number densities of 
particles and antiparticles, the condition T <C T c is most easily satisfied in the early universe. 
Let us explore in more detail this scenario for the case of DM particles. 


III. THE CONDENSATE AS A VIABLE CANDIDATE FOR DARK MATTER 

A potential CDM candidate requires, among other things, of an energy density with the 
right amount, p(a = 1) = pcDM.now, that is already decoupled from the thermal plasma at 
BBN and redshifts like a nonrclativistic component, p ~ 1/a 3 , at matter-radiation equality. 
It is well known that particles with a mass less than about an eV that are produced thermally 
in the early universe cannot contribute significantly to the present matter content. This is 
a consequence of the Cowsik-McClclland bound |56j, of which standard model neutrinos are 
a particular example. Furthermore, the energy density in these low mass thermal particles 
redshifts with the cosmological expansion as radiation, p ~ 1/a 4 , for a too large period 
exceeding the time of matter and radiation equality. They are therefore not a good prospect 
for CDM. Let us see how this picture is modified in presence of an asymmetry. 

For an asymmetric component the interactions with the thermal plasma can accommodate 
a significant amount of charge density in the ground state zero-mode. The particles in a 
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condensate behave like standard CDM, at least at the level of the background universe, and 
as long as the temperature is not close to the critical temperature of condensation so that 
the charge in the zero-mode does not change substantially with cosmic expansion, %{T) « q 
(we already argued that the condensation is a very sharp process, and as we will find next 
for any sensible candidate the temperature of decoupling is always well below the critical 
temperature of condensation). In order to see this note that, according to the identities in 


Eq. (18), the particles in the condensate contribute to the energy density and pressure of 
the gas in the form p = qnn , p — 0. The charge conservation within a comoving volume, 
q ~ 1 /a 3 , guarantees then the standard density evolution of a nonrelativistic component, 
p ~ 1/a 3 . Note that since p = qm , the amount of CDM is now fixed by an asymmetry, and 
it is not the result of the relics left dynamically at thermal decoupling (we will discuss about 
the relics left by a thermal condensate later in this section). There is no problem a priori 
with low mass candidates as long as they are not too light, m > H d , so that the description 
in this paper applies and the previous simple arguments make sense. Here H d is the value 
of the Hubble parameter at thermal decoupling. Since DM should be already decoupled at 
BBN, that imposes the constraint 

m > -Hrbn ~ 10 -14 eV (19) 


on the mass of the scalar particle. Later in Section |III 2| we will increase this number to 
10 _12 eV using current bounds on the effective number of extra relativistic species. Candi¬ 
dates lighter than this value require a more detailed analysis that we leave for another paper. 
Once the particles are decoupled from the thermal plasma the occupation numbers get frozen 
and they just redshift with the volume of the universe, n^ b {a > a d ) = n^ b (a = T d )(a d /a) 3 . 
This guarantees that the condensate does not evaporate if it did not before the time of 
decoupling, even if the temperature of the universe can grow above the critical temperature 
of condensation at later times (actually the critical temperature is not even well defined if 
the gas is not in thermal equilibrium). 

Let us put some numbers to this. We can parametrize the DM asymmetry in terms of 
its ratio to the photon number density, 


hCDM = ^ n ° w — 25.505(f2cDMh 2 ) f —^ , 

\m J 


7, now 


( 20 ) 


where q - />cDM,now/(™i 3 ), S2 C dm - Pcdm ,now/ Pc, now 5 p c = 3Mp[// 2 /(87r), and n 7 = 

2£(3)T 3 / 7r 2 . Today H now = 100/i km/(sMpc), and T now = 2.7255 ± 0.0006 K [5?], with 
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Mpi m 1.22 x 10 19 GeV the Planck mass. This assumes only one component of CDM, and 
also that all the particles in that component are in a condensed phase. We will relax these 
assumptions at the end of this section. Below we will find that for an asymmetric condensed 
DM component cosmological observations demand m < 100 eV. According to PLANCK 
^cdm^ 2 = 0.1198 ± 0.0015 at 1 a (68%CL) [58], and thus for a successful candidate we need 
a mechanism that generates a fairly large DM asymmetry, ^cdm > 1CD 2 ; compare this with 
e.g. the baryon-to-photon ratio, r] bar ~ 1CT 10 [59]. We are not going to explore this mecha¬ 
nism here, but again let us stress that it is this large value of the asymmetry what makes 
it possible that light particles behave like a nonrelativistic component in spite of being in 
thermal equilibrium in the early universe. 

As long as the scalar particles are in thermal equilibrium with the constituents of the 
standard model of particle physics, they will share a common temperature. Using the 
relation T ~ 1/a ~ (1 + z) for a universe in adiabatic expansion, we can determine this 
temperature as a function of time in terms of the temperature in the CMB photons today, 2 


T(z) = 8.617 x 10 


-5 


T 

-L rn 


K 


(1 + z) eV. 


( 21 ) 


On the other hand, from Eq. (20) and the relation q rsj \ I CL (1 + z) 3 for the charge 

conservation, we can determine the charge density necessary to match the totality of the 
CDM, also as a function of time, 


q(z ) = 8.127 x 10 11 ( — ) (0 C DMh 2 )(l + ^) d eV' 

V m 


( 22 ) 


Combining Eqs. (21) and (22) we can infer when i) the transition at T ~ m from a relativistic 
to a nonrelativistic gas of bosons occurs; ii ) the transition at q ~ m 3 from a universe with a 
high charge density to another with a low charge density takes place; and in) the transition 


2 Properly speaking, from the conservation of the entropy in a universe in adiabatic expansion, s(T)a 3 = 
s(T now ), we should write T = [g*s{T n0 w)/g*s{T)] 1 ^ 3 T now (l + z), where the function g*s(T) parameterizes 
the effective number of relativistic species contributing to the entropy density s at temperature T, see 
Eqs. (23) and (24) below for details. Then, the identity in Eq. © is only satisfied after the electron- 
positron annihilation at T ~ 0.5 MeV, once the function g*s(T) gets fixed to the current value g*s(T now ). 
However, in practice and for the model we are considering in this paper 0.3 < [<7*s(Uiow)/ g*s(T)] 1 ^ 3 < 1, 
and this factor is not going to be relevant for what we want to show here. Furthermore, this extra 
factor can only reduce the actual value of the temperature with respect to that in Eq. (21), and then it 
contributes favorably to the appearance of a condensate. 
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FIG. 1: Cosmological evolution of a gas of bosons with a non-zero charge density as a function 
of the temperature of the CMB photons T for different masses of the scalar particle m. The 
shaded region indicates the presence of a condensate, T < T c . The transition from a relativistic 


to a nonrelativistic gas occurs at T ~ m in Eq. (21), whereas the transition from a universe 
with a high-charge-density to one with a low-charge-density takes place at q ~ m 3 in Eq. ( ]22[ ) . 
Here we have considered that the condensate contains the totality of the CDM, and that it is in 
equilibrium with the thermal plasma. Note that the particles that could conform a cosmological 
DM condensate should be necessarily much lighter and weakly interacting than usually expected 
for standard WIMPs, and that for the allowed range of masses the temperature of decoupling is 


always well below the critical temperature of condensation (see Sections III 1 and III 2 for further 
details). The dotted line represent the beginning of BBN and are just for orientation. This figure 
is not applicable to candidates lighter than around 10~ 12 eV since they cannot have reached the 
regime of slow expansion, H <C m, at decoupling. 


at T = T c from a condensate to a phase with no condensate happens. All these regimes are 
illustrated in Figure [TJ 

We can look at this figure as representing, for different masses of the scalar particle, 
the phases of a cosmological gas of bosons with a charge asymmetry that matches the 
totality of the CDM, whenever the bosons are in equilibrium with the thermal plasma. 
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Note that, for a fixed value of the mass of the scalar particle, the condensate is present 
until the temperature of the universe reaches the critical temperature of condensation, when 
the condensate in principle disappears: if we move from the right to the left in Figure [l] 
following horizontal lines we eventually leave the shaded region. That means that we need 
to break thermal equilibrium before reaching that point, so that the amount of charge in 
the condensate freezes out. Furthermore, in order to do not affect the abundance of light 
elements predicted in the standard cosmological scenario, the temperature of decoupling 
should be at least as high as the temperature in the universe at the beginning of BBN, 
T ~ lOMeV. Then, we can easily appreciate from Figure [l] that the condensed phase 
is necessarily in the relativistic, high density regime, i.e. prior to the beginning of BBN 
the line that divides the relativistic and the nonrelativistic regimes is always above the 


shaded region. That implies that, according to the expressions in Eq. (18), apart from the 
condensate there will be a thermal cloud of relativistic particles. This cloud also freezes at 
the time of decoupling, producing (hot) thermal relics that will be around in the late-time 
universe and could affect the observations. This is a consequence of the thermal origin of 
the condensate. 

Let us explore in more detail the cosmological consequences of these thermal relics. Since 
the DM particles decoupled being hot, we can easily compute their temperature (which is 
only well defined as long as the particles are relativistic, and moreover does not necessarily 
coincide with the temperature of the CMB photons), and then their abundance in the present 
universe. In a universe in adiabatic expansion the comoving entropy S is conserved, and 
then S = s(T)a 3 = s(T now ). Here 


2vr 2 

>(T) = -^a.s(T)T 3 


(23) 


is the entropy per unit volume, and the function 

9 .,u)= v z jpy 

i bosons ' ' i fermions ' ' 


(24) 


counts the total number of effectively massless degrees of freedom contributing to the entropy 
density at temperature T, i.e. those that are in thermal equilibrium at a temperature 
Ti(T) rrii, or that are not in thermal equilibrium anymore but decoupled and froze out 
when they were still relativistic. Here g l is the number of spin states of the particle and 
antiparticle i, and the factor of 7/8 in the second term accounts the difference between 
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fermions and bosons. From the equation of entropy conservation, S = Ssm + Relics, we 
obtain 

ggsTO V 


T{a) = 


1/3 


~\T d 


va 

for the temperature of the CMB photons as a function of time, whereas 


T re iics(a) ={-)T d 

a 


(25) 


(26) 


for the temperature in thermal relics also as a function of time. Here T d and a d are the 
temperature and the scale factor at which DM decoupled from the primordial plasma, and 
we have used the fact that the particles in the standard model and the thermal relics do not 
interact after they decoupled, so that the comoving entropies Ssm and Relics are conserved 
independently. The two expressions above are valid only after DM decoupling, with T DM (T > 
T d ) = T before that time. Here the superscript SM in g*(T) makes reference to the fact 
that this function only takes into account the relativistic degrees of freedom in the standard 
model. For the standard model of particles physics this quantity is tabulated and can be 
found in the literature, see e.g. Fig. 3.5 in Ref. ra Note that g^^T) < 106.75, and today 
g^(Tnow) = 3.909, where photons, g 7 = 2, T 7 = T, and neutrinos, g u = 6, T v = (4/ll) 1 / 3 T, 
are not in thermal equilibrium but decoupled when they were still relativistic. This upper 
bound on gf£ (T) will be useful soon. 


Combining the two identities in Eqs. (|25|) and (|26|), we obtain 

Relics 


T 


sStm 1/3 

(Td)) 


<lfs 


(27) 


where again this expression is only valid at temperatures T < T d . Note that if after DM 
decoupling a standard model particle that is in thermal equilibrium becomes nonrelativistic 
and annihilates (decreasing the value of the effective number of relativistic species g^), its 
entropy is transferred to the other particles that are still in equilibrium, but not to the scalar 
bosons, “reheating” the CMB. That is the reason for which, in general, T re ii CS ^ T. 


The expression in Eq. (27) is valid only as long as the DM particles are relativistic. Once 
they go nonrelativistic the thermal spectrum breaks down, and in general we cannot define 
a proper temperature for the boson relics. However, we can still write 

^relics g*f(T) 


n ^ 


g!¥(T d 


(28) 
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This expression comes from the conservation of the comoving number density of thermal 
relics after decoupling, N = n re n cs (T d )/s(T d ) = n ie u cs (T)/s(T), where we have used the 
identities n rehcs (T) = n 7 (T) = 2C(3)T 3 /V 2 and g*s(T)/g* s (T d ) = gff(T)/g^f(T d ). Once 
again the number of thermal relics per photon decreases whenever a particle species self- 
annihilates and disappears heating up the CMB photons. Incidentally, for the model we are 


considering in this paper the two ratios in Eqs. (27) and (28) have not changed since the 
electron-positron annihilation at T ~ 0.5 MeV. 

Note that since we are dealing with hot thermal relics they are almost as abundant 


today as the photons in the CMB, see Eq. (28) above, and they can affect the cosmological 
observations in two different ways: if they are too “heavy”, m > 10~ 4 eV, they could comprise 
too much DM; if on the contrary they are too “light”, m < 10 6 eV, they could potentially 
spoil the success of BBN. Let us now discuss these two constraints in turn. 


1. Dark matter particles of masses m > 10 4 eV 

If the mass of the scalar particle is too large, m > 10” 4 eV, the thermal relics will be 
nonrelativistic today, m > T relicSjIlow < T now ~ 10~ 4 eV, and they will contribute to the mass 
density of the universe according to 

Q relics h 2 = °- 153 ^sM (Td) (ev) ’ ( 29 ) 

where Drelics = Prelies,now/Pc,now; Prelies = ^relics ^-relics = ‘^‘Ci.^'^relics/^ ; Pc = 3AjTpj77 /(87 t), 

and today gfs{T now ) = 3.909. Note that the lighter and more weakly interacting [i.e. the 
larger g^^Td)} the scalar particles are, the less they will contribute to the matter content. 
Even in the best case scenario in which g^(T d ) = 106.75, thermal relics heavier than 
71.546 eV will contribute more than D re iecis^ 2 ~ 0.1 to the energy density of the present 
universe, and we can then safely exclude these values for the mass of the scalar bosons; 
see Figure [T| for details. Compare this with the Lee-Weinberg bound for standard cold 
WIMPs, where m > 10 GeV [60] (particles of mass 10 eV < rn < 10 GeV that are produced 
thermally and interact mainly through the weak nuclear force would represent too high a 
DM-to-photon ratio, whereas particles of mass m < 10 eV would contribute as hot DM and 
present problems with e.g. structure formation). 


20 






2. Dark matter particles of masses m < 10 6 eV 


If on the contrary the mass of the scalar particle is too low, m < 10 6 eV, the thermal 
relics will be relativistic at the time of BBN, Xbbn ~ 10 — 0.1 MeV, and they will contribute 
to the effective number of extra neutrino species in the form 

1 


AN eS = 27.114 


kW )] 4 / 3 ’ 


where AiV e ff is defined through 


7r 


Prad = '—g*(T)T 4 = 


1 + I (3.046 + AAf„) (A) 


4/3' 


P'1 ' 


(30) 


(31) 


Here p ra d denotes the energy density in radiation, and p 7 = (7t 2 /15)T 4 that in photons. 
Note that, actually, for the whole range of masses that survives the previous constraint 
in Section III 1[ m < 100 eV, the particles are relativistic at BBN, and then the dis¬ 
cussion in this section is general. As usual N ^ = 3.046 represents the effective num¬ 
ber of neutrinos in the standard model of particle physics [ST], and T v /T = (4/11) 1 / 3 
[^relics /T = (3.909/p®|l(T d )) 1 / 3 ] is the ratio of the temperature in neutrinos [thermal relics] 
to that in photons after the electron-positron annihilation. Once again, the more weakly 
interacting the scalar particles are, the less they will contribute to the energy budget, i.e. to 
the effective number of neutrino species. However, since we are now dealing with relativistic 


relics, and contrary to what happens in Section III 1 this contribution is not sensitive to 
the actual mass of the particles. 

An analysis by Cyburt et al combining observations of the primordial 4 He mass fraction, 
light element abundances, and baryon-to-photon ratio (T^+D/H^+pcmb) results in an upper 
bound of AA^ BN < 0.804 at 68%CL to the extra neutrino species at BBN [52], and then in a 
lower bound of g™{T d ) > 13.994 to the effective number of relativistic species contributing 
to the entropy density also at BBN. Using the latest data provided by the PLANCK collab¬ 
oration in a joint CMB+BBN analysis (D+Planck TT,TE,EE+lowP), A#™ B < 0.234 [58]. 
we can improve this bound to g™(T d ) > 35.317, now at 95%CL. Note however that this 
latter constraint applies only for DM bosons of mass less than about an eV, which are still 
relativistic at recombination. These values of g^g 1 lie around the scale of the QCD phase 
transition, Aqcd ~ 200 MeV, where the baryons and mesons are formed and the number 
of massless degrees of freedom drops drastically from gfg 1 ~ 60 to gfg 1 ~ 20. We can then 
conclude that T d > Aq CD for an asymmetric condensed DM candidate. 
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Two comments are in order here. First, since the number of massless degrees of freedom 
in the standard model is constrained to < 106.75, then there is a lower bound for 

the number of extra neutrino species: AA r e ff > 0.054. This value is around the sensitivity 
of ±0.04 on N e g expected from future measurements of the polarization of the CMB [63], 
and then the cosmological observations could test the viability of this model soon. Second, 
the lower bound of Td > 200 MeV on the temperature of freeze out is related to the self¬ 


annihilation, DM±DM=f±SM±SM, and elastic scattering, DM±SM^DM±SM, cross sec¬ 
tions associated to the processes that maintained the condensate in chemical and kinetic 
equilibrium with the thermal plasma in the early universe, respectively, and then to the 
strength of the interactions between DM and standard model particles. As a consequence 
of this high value of the temperature of DM decoupling the interaction of the bosons with 
the standard model particles should be very weak, even weaker than in the case of ordinary 
WIMPs, where e.g. Td ~ 1 MeV for neutrinos. 3 Note also in Figure [l] that for the allowed 
range of masses of the scalar particle the critical temperature of condensation is always above 
the temperature at which thermal equilibrium broke down, and also above the temperature 
of BBN. This would guarantee the survival of the coherent state to thermal evaporation, 
however, this point deserves a more detailed analysis of the kinetic and chemical decoupling 
in terms of the Boltzmann equation in the presence of a condensate. 

Let us come back to the general picture. One could imagine scenarios in which the conden¬ 
sate would represent only a fraction of the total CDM in the universe, the remainder being 
in the form of standard thermal relics or any other candidate to DM. In this scenario the 
lines in Figure [l] dividing the high-density/low-density, and the condensed/non-condensed, 
phases move down (note that they are sensitive to the asymmetry of the universe, and in 
particular they disappear if q — 0), but the other information in the figure remains unaltered. 
This does not affect, in particular, the conclusion that the thermal relics accompanying the 


condensate are hot, nor the upper bound we found previously in Section III 1 for the mass 
of the scalar particles. Thermal relics of mass m < 100 eV would contribute to the energy 


density of the universe as hot DM, and they cannot represent a significant part of the matter 


3 For the case of standard cold thermal relics the temperature of decoupling can be also as high as T d ~ 
lOOMeV, but there this is a consequence of the rarefaction of the gas due to the unbalanced self-annihilation 
of DM into standard model particles, and not of the weakness of the interactions. 
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content. 4 Then, if there are no more particles in the universe apart from the DM boson and 
the constituents of the standard model, we can conclude that the totality of the CDM should 
be necessarily in a light condensate, or in the form of heavy thermal relics, but we cannot 
have a combination of both. This last statement could be avoided, of course, if we had 
another component of DM (apart from the boson particles that conform the condensate) in 
the form of e.g. WIMPs, but the previous constraints on the scalar particles still apply. 

To summarize, the scalar DM scenario we presented in this paper has necessarily a con¬ 
densed and a thermal component. The condensed component is essentially determined by 
a charge asymmetry, which must be substantial to provide enough CDM. Lest the thermal 
component represent too much DM (or the right amount, but with wiping out too much 
small-scale structure), it must be lighter than about a hundred of eV’s, and then the DM 
boson must decouple at high temperature so that it does not add a light degree of free¬ 
dom during BBN. These constraints remain even if the scalar DM particles comprise just a 
significant part, rather than all, of the observed CDM. 

* * * 

Since early decoupling is crucial, it is worth some discussion here. The question of 
decoupling is intimately related with the question of what led the system into thermal 
equilibrium. For a standard thermal candidate it is usually argued that gravity is not 
sufficient, and there should be an additional interaction, usually the weak nuclear force, that 
couples DM to the particles in the standard model. Here the argument goes as follows [T]: 
For gravity the interaction rate scales like r grav ~ np>yi(<jv) ~ T 5 /Mp 1 , where «dm ~ T 3 
and (cry) ~ T 2 /Mpj. On the other hand, during the radiation dominated era the Hubble 
parameter is proportional to H ~ T 2 /Mpi, and then decoupling, r grav ~ H , occurs at 
Thee ~ Mp\. We cannot trust these expressions at the Planck scale, and then we conclude that 
the gravitational interaction alone could not have brought the early universe into thermal 
equilibrium. 

However, for the case of an asymmetric DM component the number of particles scales with 

4 Hot thermal relics free-stream in the early universe. This avoids gravitational instability and induces a 
cutoff in the mass power spectrum at the scale Mj ~ Mpim~ 2 ~ (10keV/m) 2 x 1O 8 M 0 |M]. In order to 
not erase dwarf galaxies, Mdwarfs ~ 1O 9 M 0 , and heavier structures we should demand m > lOkeV. 
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tt-dm ~ T 3 + q ~ (1 + ? 7 cdm)F 3 , and we can then conclude that T d ~ (1 + ?7cDM) _1 ^ 3 d^pi. For 
ultralight candidates, ^cdm 3> 1, the interaction rate increases with respect to the standard 
scenario and decoupling is expected to take place well after the Planck era, where we can 
now rely on these expressions. This is of course only a very naive picture; a more accurate 
treatment would employ the Boltzmann equation [63] to treat the decoupling process in the 
presence of a condensate. We leave this analysis for a future work. 


IV. CONCLUSIONS 

In this paper we have considered the possibility that much or all of the CDM in the uni¬ 
verse results from an asymmetry between the number densities of particles and antiparticles, 
rather by the abundance of thermal relics frozen out at the time of thermal decoupling. We 
have identified three conditions that, if satisfied, guarantee the appearance of a thermal 
BEC during the early stages of the universe: i) DM is described in terms of a boson field 
with a conserved current; ii ) There is a nonvanishing total boson charge in the universe; and 
in) The boson was in thermal (kinetic and chemical) equilibrium with the particles in the 
standard model. If, in addition, iv ) Thermal decoupling took place at an energy scale below 
the critical temperature of condensation, then the macroscopic excitation of the ground state 
freezes out and could survive the ensuing cosmological evolution until now. Note that for 
the purposes of this paper we did not need to develop any boson candidate, not even specify 
whether they are fundamental, or just an effective low energy description, e.g. made of pairs 
of fermions. We have estimated the charge asymmetry required to account for the totality 
of the CDM using the condensate, but do not address here how this large asymmetry would 
have arisen. 

According to our findings there is an upper limit on the mass of the DM particles in a 
thermal BEC. This is because in this model the macroscopic excitation of the zero-mode 
has a relativistic origin and then, apart from the condensate, there will be thermal relics 
in an abundance close to that of the photons in the CMB. These particles must then have 
masses lower than around 100 eV so that their density is below the current cosmological 
constraints (compare this with the unitarity bound m < 100 TeV one finds in the case 
of ordinary WIMPs |66j). The existence of thermal relics could also alter the success of 
BBN. In order to not contribute in excess to the effective number of neutrino species the 
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relies should be cold enough, that is, they should decouple early on in the history of the 
universe. Using the latest BBN and BBN+CMB observational data we obtain that the DM 
bosons should be already decoupled at the time of the QCD phase transition, when the 
temperature of the universe was of the order of 200 MeV. As a consequence the interaction 
between the DM particles that conform the condensate and those in the standard model 
should be even weaker than the weak nuclear force. We speculate about the possibility that 
gravity alone could be the responsible of the early thermal equilibrium. Incidentally, there 
is a lower bound to the excess of neutrino species predicted by an asymmetric condensed 
DM component, A N e g > 0.054. This value could be resolved in the next future using CMB 
data, and then test the viability of the model. 

Let us emphasize that, in spite of the early thermal equilibrium, the particles in a DM 
condensate could be as light as m ~ 10 -12 eV and behave like a nonrelativistic component 
in terms of their background “equation of state.” (For candidates of masses lower than 
this value the decoupling from the thermal plasma necessarily occurs during the regime 
of fast expansion, H m, and then a more elaborate analysis is necessary to treat the 
phenomenology of ultralight particles.) We usually tend to think of thermal cold candidates 
as particles with large masses, as happens in e.g. the standard symmetric cosmological 
scenario. However, if DM is described in terms of a boson field, there is another possibility 
to cool the particles in the universe: giving them a large asymmetry at the expense of low 
masses. 

Note that even though the input of the asymmetric condensed DM scenario is very sim¬ 
ilar to that of ordinary WIMPs (i.e. a thermal quantum held theory), the output differs 
drastically: whereas at the effective level we can think of WIMPs as a collection of classical 
particles, the condensate is described in terms of a classical field theory. This can affect, for 
instance, the process of structure formation at small scales, leading to possible interesting 
observational consequences in the late-time universe [SUSIEM. 

The analysis in this paper complements the standard thermal scenario of a boson DM 
candidate (without an asymmetry, or with an asymmetry much lower than the value reported 


in Eq. (20), so that the quantum Bose statistics is blurred by the classical Maxwellian 
distribution for most of the history of the universe), and indicates that a boson particle that 
was in thermal equilibrium with the constituents of the standard model of particle physics 
should be in a condensate, or in the form of thermal relics, but we cannot have a combination 
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of both aspects contributing significantly to the matter content today. 
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Appendix A: Some properties of the polylogarithm functions 


In this appendix we present some useful limiting values of the polylogarithm functions. 


They are useful when evaluating the expressions in Eqs. (16) in the main text. 

.2 


7T 


Li 3 (x —> 1) — C(3) + —(x — 1) + 
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7r 


Li 4 (a; -A 1) — — + C(3)(o; - 1) + 


(Ala) 
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Here ((x ) is the Riemann zeta function, with £(3) ~ 1.2021. 
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